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The frequently discussed strongly interacting gauge theory with a fermion flavor doublet in the 
two-index symmetric (sextet) representation of the SU(3) color gauge group is investigated [1]. 
The chiral condensate and the mass spectrum are consistent with chiral symmetry breaking (^SB) 
at vanishing fermion mass. In contrast, sextet fermion mass deformations of spectral properties 
are not consistent with leading conformal scaling behavior near the critical surface of a conformal 
theory. A recent paper [2] which could not resolve the conformal fixed point of the gauge coupling 
from the slowly walking scenario of a very small nearly vanishing j3 -function is not in conflict 
with jSB reported here. A light Higgs impostor could emerge as the dilaton from spontaneous 
symmetry breaking of scale invariance or, without the dilaton mechanism, as a composite state. 
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1. Introduction 

The new Higgs-like particle with decay modes not far from that of the Standard Model brings 
new focus and clarity to the search for theoretical frameworks. One example is the light dilaton as 
a pseudo-Goldstone particle of spontaneous breaking of scale invariance that has been featured in 
recent phenomenological discussions as a viable interpretation of the discovery. Nearly conformal 
gauge theories serve as theoretical laboratories for credible realizations of this scenario (this short 
report is based on [1] which provides references to the history of the field including recent work). 
We investigate here a candidate theory with a fermion flavor doublet in the two-index symmetric 
(sextet) representation of the SU(3) color gauge group close to the conformal window with a small 
beta function, as suggested by results of a recent paper [2]. With the limited statistical accuracy of 
the nearly vanishing /3 -function the existence of a conformal fixed point gauge coupling remains 
unresolved from the alternative slowly walking scenario. Consistency with #SB, reported here 
for small fermion mass deformations, would require the sextet model to remain just below the 
conformal window with a very small non-vanishing /3-function (see, also [3,4]). In this case 
the model which exhibits the simplest composite Higgs mechanism leaves open the possibility of 
a light scalar state emerging as the pseudo-Goldstone dilaton state from spontaneous symmetry 
breaking of scale invariance. Even if scale symmetry breaking is entangled with #SB without 
dilaton interpretation, a light Higgs-like composite state can emerge close to the conformal window. 
We outline a new lattice Higgs project to resolve these important problems. 

2. Computational strategy and lattice simulations 

Probing #SB, and conformal behavior for comparison, we extrapolate the spectrum to infinite 
volume at fixed fermion mass m. In large volumes the leading finite size corrections are expo- 
nentially small and dominated by the lowest state of the spectrum which is expected to have pion 
quantum numbers. From the mass spectrum, extrapolated to infinite volume, we can probe the pat- 
tern of #SB when small fermion mass deformations are simulated close to the massless limit. We 
also probe the hypothesis of mass deformed conformal scaling behavior and find results strongly 
favoring the #SB hypothesis. 

A new analysis is presented at gauge coupling jS = 3.2 based on a subset of runs in the fermion 
mass range m = 0.003-0.010 on 24 3 x 48, 28 3 x 56, and 32 3 x 64 lattices. Five fermion masses 
at m=0.003, 0.004, 0.005, 0.006, 0.008 are used in most fits probing the two hypotheses. A very 
large and expensive 48 3 x 96 run was added recently at m = 0.003 to control finite size effects. 
Simulation results at jS = 3.25 were also obtained in the mass range m = 0.004 — 0.008 on 24 3 x 48, 
28 3 x 56, and 32 3 x 64 lattices. We have used the tree-level Symanzik-improved gauge action with 
stout smeared staggered fermions for all results reported here. The normalization of the lattice 
coupling j3 and details of the simulations are discussed in [1]. 

To control finite size effects, infinite-volume extrapolations were performed for the lowest 
state in the spectrum with pion quantum numbers, the related decay constant F % , and the chiral 
condensate (\j7y) ■ They are shown in Figure 1 where g\ (X , 7] ) describes finite volume corrections 
from the exchange of the lightest pion state with X = M K L and lattice aspect ratio r\ = T/L, sim- 
ilarly to what was introduced in [5]. The fitting procedure approximates the leading treatment of 
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M p = 3.2 m=0.003 F p = 3.2 m=0.003 (3>B) m=0.003 finite volume fit 
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Figure 1: Finite volume dependence at the lowest fermion mass for j8 = 3.2. The form of g\{X,r\) is a complicated 
infinite sum which contains Bessel functions and requires numerical evaluation [7], Since we are not in the chiral log 
regime, the prefactor of the gi(A,Tj) function was replaced by a fitted coefficient. The leading term of the function 
gi (A, r]) is a special exponential Bessel function K\ (A) which dominates in the simulation range. 

the pion which wraps around the finite volume, whether in chiral perturbation theory (#PT), or in 
Liischer's non-perturbative finite size analysis [6]. This equivalence relaxes the requirement on the 
fitted parameters cm,cf,c\ to agree with 1-loop #PT as long as the pion is the lightest state domi- 
nating the finite volume corrections. The infinite-volume limits of M K , F n , and for m = 0.003 
at j6 =3.2 were determined self-consistently from the fitting procedure. Similar fits were applied 
to other composite states. The value of M n in the fit of the left plot in Figure 1 was determined 
from the highly non- linear fitting function and used as input in the other two fits. Based on the 
fits at m = 0.003, the results are within one percent of the infinite-volume limit at M n L = 5. In the 
fermion mass range m > 0.004 the condition M n L > 5 is reached at L = 32 and we will interpret the 
results from the 32 3 x 64 runs for m > 0.004 as infinite-volume behavior in mass deformed chiral 
and conformal analyses. 



3. The chiral condensate 



Our simulations show that the chiral condensate (yy) is consistent with ^SB and remains 
non-vanishing in the massless fermion limit. The cutoff dependent UV contributions are identified 
at finite fermion mass. The linear mass term c\ (a) ■ m is a quadratically divergent UV contribution 
w aT 2 • m with lattice cutoff a. There is also a very small third-order UV term cs(a) • m without 
power divergences which is hard to detect for small m and has not been tested within the accuracy 
of the simulations. IR finite contributions are identified in the chiral expansion of the condensate. 
There is an m- independent constant term which is proportional to BF 2 , a linear term proportional to 
B 2 ■ m, a quadratic term ~ B 3 F 2 ■ m 2 , and higher order terms, in addition to logarithmic corrections 
generated from chiral loops. The expansion in the fermion mass is expressed in terms of low energy 
constants of chiral perturbation theory, like B and F [8]. A more complete description of our fitting 
functions is given in [1]. We used two independent methods for the determination of the chiral 
condensate in the massless fermion limit. In the first method fits were made directly to with 
constant and linear terms in the fitted function. Quadratic and third order terms are hard to detect 
within the accuracy of the data. The result is shown in the left plot of Figure 2. When the quadratic 
term is added to the fit, the massless intercept Co = {Yy) m =o from the quadratic fit agrees with 
the one from the linear fit and the quadratic fit coefficient in C2 ■ m 2 is zero within fitting error. For 
an independent determination, we also studied the subtracted chiral condensate operator defined 



3 



The sextet gauge model, light Higgs, and the dilaton 



Julius Kuti 




Figure 2: The chiral condensate and its reduced form with subtracted derivative (both have to converge to the same 
chiral limit) are shown in the left plot with linear fit to the condensate. The data without derivative subtraction cannot 
detect higher order fermion mass terms with significant accuracy. The fit to the reduced form with subtracted derivative is 
defined in [1] and shown in the magnified right plot. A linear term is not included in this fit since the subtracted derivative 
form approximately eliminates it. The value of do at in = is shown to be consistent with the direct determination of 
Co from the chiral limit of (xjfyf) . The consistency is very reassuring since the two results are derived from independent 
determinations. For m = 0.003 the data from infinite-volume extrapolation were used in the fit. As we explained earlier, 
at higher m values the largest volume 32 3 x 64 runs were used for the condensate and its derivative subtraction. 



with the help of the connected part % C onn of the chiral susceptibility %. Once the derivative term 
is subtracted, the first non-perturbative IR contribution, quadratic in m, is better exposed. The two 
independent determinations give consistent non- vanishing fit results in the massless chiral limit as 
shown in Figure 2. 

4. Spectral tests of the #SB hypotheses 

The chiral Lagrangian describes the low energy theory of Goldstone pions and non-Goldstone 
pions in the staggered lattice fermion formulation. It can be used as an effective tool probing the 
%SB hypothesis at finite fermion masses including extrapolation to the massless chiral limit. Close 



sextet modGl Goldstone pion in PCAC channel sextet model from PCAC channel 




Figure 3: Polynomial fits from the analytic mass dependence of the chiral Lagrangian without logarithmic loop correc- 
tions are shown for the Goldstone pion and F n . The dashed line in the left plot for the Goldstone pion shows the leading 
linear contribution. The data point at m = 0.003 is determined from the 48 3 x 96 lattice and the other fitted points are 
based on the 32 3 x 64 runs. 

to the chiral limit, the pion spectrum and the pion decay constant F % are organized in powers of 
the fermion mass m which is an input parameter in the simulations. Chiral log corrections to the 
polynomial terms are generated from pion loops [9]. Their analysis will require an extended dataset 
with high statistics. 
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In Section 2 we presented results of infinite-volume extrapolations. Based on these observa- 
tions, in fits to the observed pion spectrum and F„, we will use infinite-volume extrapolation at 
m = 0.003 and treat the 32 3 x 64 runs for m > 0.004 as if the volume were infinite. In Figure 3 we 
used the local pion correlator with noisy sources to extract M n and F n . The correlator is tagged as 
the PCAC channel since the PCAC relation, based on axial Ward identities, holds for this correlator 
and the decay constant F K can be directly determined from the residue of the pion pole. Based on 
the analytic fermion mass dependence of the chiral Lagrangian, and using the lowest four fermion 
masses, good polynomial fits were obtained for M % and F % as shown in Figure 3 with fitting func- 
tions M 2 . = IB ■ m + pm ■ m 2 and F n = F + pp ■ m . The parameters B and F are defined in the 
two leading terms of the chiral Lagrangian [9]. In this simple fitting procedure B is set from the 
fit to M 2 and F is set independently from the fit to F % . The fit parameters pm and pp describe 
the respective leading power corrections in m to the pion mass and decay constant. Limited to a 
single lattice spacing only, the accuracy of our dataset is not sufficient for the robust determination 
of chiral log corrections that will provide important consistency check for B and F in our future 
analysis of an extended dataset. This analysis requires rooted and partially quenched staggered 
perturbation theory at finite lattice spacing for simultaneous fits of M n and F % with a consistent 
pair of cutoff-dependent F and B values [10]. We made the first step in this direction by adding a 
new run set to our database at jS = 3.25. Reduction in taste breaking is significant at /3 = 3.25 with 
smaller lattice spacing. Our staggered perturbation theory analysis will be presented elsewhere. 

5. Spectral tests of the conformal scaling hypothesis 

Under the conformal scaling hypothesis, the mass M n and the decay constant F % are given at 
leading order by M n = cm ■ m 1 ' 1+y and F % = cp ■ m l l l+y . The coefficients cm and cp are channel 
specific but the exponent y is universal in all channels [11]. The leading scaling form sets in for 
small m values, close to the critical surface. According to the hypothesis, there is an infrared 
conformal fixed point on the critical surface which controls the conformal scaling properties of 
small mass deformations. All masses of the spectrum can be subjected to similar conformal scaling 
tests, but we will mostly focus on accurate data in the M K and F % channels. 

When M n and F % are fitted separately in the range of the four lowest fermion masses closest 
to the critical surface, we get reasonable % 2 values for the fits. However, the incompatibility of 
the fitted y values disfavors the hypothesis, inconsistent with mass deformed conformal behavior. 
The conflicting simultaneous fits to universal conformal form with the same y for the Goldstone 




Figure 4: The first plot shows the simultaneous conformal fit result for the pion mass, while the second displays the 
M K residuals. The last two plots show the simultaneous fit result for the pion decay constant and the F n residuals. The 
combined fit forces y = 1.53(28) with an unacceptable ^ 2 /dof of 44.5. 
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pion and the F % decay constant are illustrated in Figure 4. Fitting to the pion mass separately 
requires 7 = 1.040(73) while the separate F % fit is forcing 7 = 2.20(15). In the combined fit they 
compromise with 7 = 1.53(28) and the unacceptable # 2 /dof of 44.5. It is important to note that 
the exponent 7 for the fit to M n only is what #SB would prefer. The separate conformal exponent 
7 for F n is large to force to the origin the linear string of data which extrapolate to a finite constant 
in ^SB. This creates conflict with the universal exponent 7 in the conformal analysis. 

From the tests we were able to perform, the sextet model is consistent with ^SB in leading 
approximation and inconsistent with conformal symmetry. It will require further investigations to 
show that subleading effects cannot alter this conclusion. It is particularly important to strengthen 
the analysis of the #SB hypothesis reaching the necessary accuracy to identify chiral log effects. 
We will also consider comprehensive conformal finite size scaling (FSS) tests which do not rely 
on infinite-volume extrapolation in the scaling fits. It remains difficult to reconcile #SB and large 
exponents in the fermion mass dependence with the low value of 7 defined by the chiral condensate 
using the Schodinger functional for massless fermions [2]. 



6. The new sextet Higgs project 

Figure 5 shows the fermion mass dependence of the scalar ++ (/o) meson without including 
the disconnected part of correlator I in Table 1 of [12]. The non-Goldstone scPion and fo are 
opposite parity states in this staggered correlator. The quantum numbers of the fo meson match 
those of the ++ state in the staggered correlator. Close to the conformal window the scalar ++ 
meson is not expected to be similar to the a particle of QCD. If it turns out to be light, the full scalar 
state, including the disconnected gluon annihilation diagram and mixing with the ++ glueball 
state, could replace the role of the elementary Higgs and act as the Higgs impostor. It is very 
difficult to do the full calculation including the disconnected gluon annihilation diagram and mixing 
with the ++ glueball state. This is the main part of our next generation sextet Higgs project. 



effective mass Ivf : ° rl from ++ connected correlator 




0.005 0.01 0.015 

m 



Figure 5: The linear fit is shown to the mass of the scalar ++ (/o) meson from the connected part of correlator I in 
Table 1 of [12]. For comparison, the scPion which is the opposite parity state in the staggered correlator is also plotted 
(at m = the mass of the non-Goldstone scPion will vanish in the continuum). Gluon annihilation and mixing of the 
fo state with a low mass ++ glueball state could generate a light scalar mass close to the conformal window. The 
disconnected part of the staggered correlator and its mixing with the glueball correlator are required to resolve this issue. 

The linear fit from the connected diagram is shown in Figure 5. It has a non-zero intercept in 
the chiral limit with a mass more than five times F so it corresponds to a heavy state and not a light 
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Higgs candidate. Close to the conformal window it is reasonable to expect that the disconnected 
diagram, after mixing with the ++ glueball state will lead to a light scalar mass in F units when 
the chiral limit is taken. The light scalar state after mixing remains a viable Higgs candidate. Only 
new simulations can resolve the issue and perhaps eliminate this attractive scenario. 

To study flavor-singlet mesons, we need to consider fermion loops which are disconnected (of- 
ten called hairpin diagrams). Flavor-singlet correlators have fermion-line connected and fermion- 
line disconnected contributions from the hairpin diagrams. To evaluate disconnected quark loops 
with zero momentum, we need to sum over propagators from sources at each spatial location for 
a given time slice. To avoid the very costly &(V) inversions to compute all-to-all propagators in 
lattice terminology, random sources have to be used with noise reduction. 

A very interesting remaining challenge and complication is the existence of two types of dis- 
tinct ++ scalar states. One of them is the composite meson state and the other is the scalar glueball 
with the same ++ quantum numbers. In dynamical sextet simulations, these two types of state will 
mix producing two scalar states. To resolve this will require a well-chosen variational operator set 
with room left for a light scalar state to emerge in the spectrum. It is also entirely possible that 
careful lattice calculations will shut down the light Higgs interpretation. 
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